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ABSTRACT 

One of the crucial aspects of density perturbations that are produced 
by the standard inflation scenario is that they are Gaussian where seeds 
produced by topological defects tend to be non-Gaussian. The three point 
correlation function of the temperature anisotropy of the cosmic microwave 
background radiation (CBR) provides a sensitive test of this aspect of the 
primordial density field. In this paper, this function is calculated in the 
general context of various allowed non-Gaussian models. It is shown that by 
COBE and the forthcoming South Pole and Balloon CBR anisotropy data 
may be able to test provide a crucial test of Gaussianity. 



PACS number: 98.80.Bp, 98.80.Dr 



Testing for the Gaussianity of the primordial fluctuation spectrum is of 
critical importance to many cosmological models. In particular, traditional 
cosmic inflation specifically predicts a Gaussian density fluctuation spec- 
trum. The scale invariant quantum fluctuations generated during the infla- 
tionary epoch are expected to serve as the primordial density perturbations 
which develop into the large scale structures we observe today [§. Compet- 
ing models for structure formation, including topological defects originating 
from cosmological phase transitions and non-standard inflation models]^], 
will also generate a scale invariant (or nearly scale invariant) power spectrum 
for density perturbations similar to that of inflation. However, the statistics 
of these latter fluctuations are non-Gaussian. Thus, the Gaussianity of the 
fluctuations provides a unique handle in discriminating different structure 
formation scenarios. In this letter, we will discuss how to test this aspect 
of the primordial density field through the temperature anisotropy of the 
cosmic microwave background radiation (GBR). 

As we showed [0, in momentum space, the lowest order deviation from 
Gaussianity is described by the bispectrum of the gravitational potential 0, 
P(t>iki,k2) =< (f>ki<Pk2'P-ki~k2 >■ When the perturbation is adiabatic so that 
the temperature anisotropy is related to the gravitational potential at the 
last scattering surface through the Sachs- Wolfe formula: 

T 3' ^ ' 

the three point temperature correlation function is related to the bispectrum 
through 



where rjo = 2Hq is the distance to the last scattering surface, and m, n, / 
are the beam directions. A non-vanishing three point function clearly indi- 
cates that the bispectrum is not zero. Note that for Gaussian primordial 
perturbations, the bispectrum is strictly zero in all cases. Thus, the three 
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point temperature correlation function is a clean test of the Gaussianity of 
primordial fluctuations. 

In 0, Falk et al. found that in an inflationary model with cubic self- 
interaction, P(f,{ki, k2, ks) is given by: 



where f3 ~ 10^^. In this paper, we will show that without invoking any new 
assumptions about inflationary models, the non-linear gravitational evolution 
of the initial Guassian perturbations will give rise to a three point correla- 
tion function, which has a similar angular dependence to that which certain 
non-Gaussian inflationary models predict, but a much larger amplitude than 
the one considered in [0. Then, we extend the analysis in to include 
more general cases of inflation, which produce not only the scale invariant 
but also the "tilted" perturbation spectrum. The extended analysis is helpful 
in discussing the effect of spectral index on the angular dependence of the 
three point function. To choose different non-Gaussian inflationary models 
through three point temperature correlation function will be hard because of 
the gravitational evolutionary effects. However, the three point correlation 
function produced by a cosmological phase transition tends to have a dis- 
tinctive angular dependence, which should enable one to prove or disprove 
the scenario through observations. Finally, we briefly discuss the effect of 
noise in the sky signal of GBR measurements on the analysis of three point 
temperature correlation function. 

By taking into account non-linear gravitational evolution, it is found that 
there are two terms which contribute to the GBR temperature anisotropy: 



where the flrst term is the Sachs- Wolfe @] effect due to the gravitational 
potential at the last scattering surface, and the second term is the generalized 
Rees-Sciama effect P due the evolution of the gravitational potential along 
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the photon path. When we adopt a flat cosmological model {Q = 1), the 
quasi- nonlinear analysis gives p|, |10|, ^: 



<P{k, 7]) = <P,{k) + a{r]) J J{k, k', k - k')<Pi{k )<j)i{k -k), (5) 

where (pi is the gravitational potential at the last scattering surface, 0(77) = 
is the expansion factor of the universe after decoupling and 

Jik, I, m) = 2{l ■ m) + ^ + \^ ■ (6) 

We first estimate the amplitude of the second term relative to the first term: 
since the expansion factor a after decoupling is ~ {1 + Zd ec) ~ 1000, the am- 
plitude of the gravitational potential at the last scattering surface is around 
10~^ as suggest by COBE, thus the ratio of the Rees-Sciama term to the 
Sachs- Wolfe term is of order 0.01-0.1. As the non-linear effects are con- 
tained in the Rees-Sciama term, it is this term that contributes significantly 
to the three point correlation function. For comparison, the non-linear term 
considered in is 10~^ times smaller than the linear term. This is a poten- 
tial problem on testing inflationary models through three point temperature 
correlation function. To be observable the amplitude of the non-Gaussianity 
produced in these models has to be large enough so that the gravitational 
evolution cannot completely dominate. 

The generic form of the bispectrum in inflationary models is given by: 

P4>{K h, ks) = \[P^{ki)P^{k2) + P^ik2)P^{k3) + P^{k3)P^iki)], (7) 

where A is a constant and =< (f){k)(j){—k) > is related to the power 
spectrum P{k) simply through P(f,{k) = P{k)k~'^. The cubic self-interaction 
model corresponds to the case where A ~ 10^^ with a scale invariant density 
perturbation spectrum, or equivalently, P(f,{k) ~ k~^; the non-linear gravita- 
tional evolution effect corresponds to the case where A ~ 2(1 + Zdec)/^ ~ 200 
and P<i,{k) ^ k'^. 
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Note that the two point temperature correlation function is related to P^p 
through: 

The three point correlation function calculated from the bispectrum given 
by Eq.^ is: 

^T{m, h, i) = 3A[C2(m, h)C2{n, i)+C2{m, n)C2{m, /) +C2(m, l)C2{h, I)]. (9) 

This is a theoretical relation between three point function and two point 
function since the finite beam size effects haven't taken into account yet. 
The formal treatment of the finite beam effect in CBR experiment can be 



found in [jT2|, |T3|. The beam can be well approximated as a Gaussian: 

/(|m-n|,a) = ^e-l™-'^l'/2<^^ (10) 

and the observed temperature correlation function will be the convolution of 
the theoretical correlation (infinite thin beam) with the beam, which is 

Cs{\m,h,i\,a) = 

(ifij^fifig'^^s/d'^ ~ ""^'1) '^)f{\^ ~ ''^'l) '^)f{\^~ ^1) '^)C3{\m' , \n' , k'\, 0). (11) 

For a special configuration of three beams where m-n = h- l = l- m = 
cos a, the beam-smoothed three point function is well approximated^ as 
[C2(cosa|cr)]^ where [C2(cos a|cr) is the two point function with the monopole, 
dipole and quadruple terms removed. Since the three point function is the 
products of two 2-point functions, it has a stronger dependence on the power 
spectra index n. Mutipole expansion of the 2-point function gives: 

C2(m, n) = J2 Ci{2l + l)Pi{m ■ h). (12) 
I 



For a power law spectrum P{k) ~ fc", Ci is given by 

^ _ i, g,^., ,, r(2/ + n-i) r(gf^, 

' 5^ To ^ r(2/ + 5-n) ■r(^^' ^ ' 
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where Qrms is the COBE measured quadruple |TT[] and To is the black-body 
temperature of CBR. From the analysis of two point correlation function, 
COBE can only put a loose bound on power spectra index |[Tl[]: n = 1.1 ±0.5. 
In Fig. 1, we plot the three point function for two different power spectra: a 
scale invariant n = 1 spectrum and a "tilted" spectrum where n = 0.7. Notice 
how the three point function depends strongly upon the power spectra index. 
Thus it is anticipated that the analysis of the three point correlation function 
will put a more stringent bound on n. 

In order to test cosmological structure formation scenarios through the 
three point temperature correlation function, we should have a clear handle 
on what various models predict. In the following, we will show that the 
cosmological phase transition can produce distinctive angular dependences 
other than the form given above. 

Cosmological phase transitions are widely discussed in the context of the 
structure formation |14]. In the case of a primordial phase transition, the 
horizon size at the epoch of phase transition is small and topological defects 
will form according to the Kibble mechanism |15| . The analysis of the three 
point correlation for defect-induced temperature anisotropy depends crucially 
on the evolution of the defect-network and the work along this line is still 
in progress. In this paper, we will show that it is instructive to consider 
initially the three point function in the late-time phase transition (LTPT) 
scenario |T6[. The calculation is considerably simplified in LTPT models 
for the following reasons: (1) the last scattering surface is assumed smooth 
in LTPT models, thus temperature anisotropies are solely produced by the 
generalized Rees-Sciama effect, 

since the fluctuations in density and gravitational potential are generated by 
the critical fluctuations at the critical point of the phase transition, ^ = 
4>S{rj — r]p), where r]p is the conformal time at the phase transition point. 
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Thus, in this latetime phase transition model, the temperature anisotropy 
takes the following simple form: 

^ = 20,=,,. (15) 

(2) For LTPT, the horizon size is large so that the finite horizon-size effect 
is negligible. We can calculate the three point correlation function from 
symmetry considerations. As pointed out by Polyakov [|l7l, the three point 
correlation function of the fluctuating field is completely determined up 
to a dimensionless constant by the conformal symmetry of the system at the 
critical point. The explicit form for the three point function is given by 

6 =< 'ip{xi)ij{x2)ij{x3) >= r]C2{xi, X2)C2(X2, X3)c2{x3, Xi), (16) 

where C2{xi,X2) =< ip{xi)'ilj{x2) > is the two point function and r] is a 
constant. In this letter, we assume that the gravitational potential is 
directly proportional to the underlying fluctuating field ip. For this case, the 
three point temperature correlation function has the following simple relation 
to the two point function: 

^T{m, n, l)=A- C2{m, n)C2(n, /)C2(m, /), (17) 

where A is a dimensionless constant. The full beam-smearing effects of the 
three point correlation function given above is messy and we will report 
it elsewhere. However, in the special case when rh-h = n- l = l- m = 
cos a, it can be approximated as [C2(cosQ;|cr)]'^, where C2(cosa|cr) is the 2- 
pt function with finite beam width a, with monopole, dipole and quadruple 
terms subtracted. We plot the approximated three point function generated 
by the phase transition in Fig. (2). 

The result obtained from Eq. (9) & Eq. (13) strongly suggest that the 
general form of the three point function, expressed in terms of two point 
functions, is given by: 

iriTTi, n, i) = Q ■ [C2(m, h)C2{n, 1) + C2(m, n)C2(m, /) + C2(m, i)C2{n, I)] 

+A ■ C2{m, n)C2{h, /)C2(m, f),(18) 
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where Q and A are constants. This is the archetype form of the three point 
correlation function that the experimental analysis should be compared with. 

The Gaussianity can be tested through the existing COBE and the forth- 
coming South Pole and Balloon GBR anisotropy data by three point temper- 
ature correlation function. In this letter, we have focused on the GOBE data 
although the idea and method discussed can equally apply to South Pole and 
Balloon experiments. The data set from COBE DMR is especially suitable 
for carrying out this test. On the one hand, the beam width of COBE is 7°, 
which is much larger than the horizon size at decoupling (~ 2°). Most non- 
linear causal processes which may lead to non-Gaussian signatures on the 
cosmic microwave background (GMB) sky are smoothed out by the beam. 
On the other hand, the COBE GMB map covers the whole sky. Thus, the 
boundary effects will be minimized. However, the detected sky signal contain 
both the intrinsic GBR temperature fluctuation and the instrumental noises, 

5T _6T ST 

TjT/obs TjT/CBR ~l 'JjTJnoise- v J 

The signal to noise ratio of the GOBE data is 1:1 and this is typical in all 
current GBR temperature anisotropy experiments. Thus, it is important to 
consider the noise term seriously in the analysis of the three point correlation 
function. Even if future analysis of the GOBE data do find a non-vanishing 
three-point temperature correlation, it may due to the instrumental noise. 
However, if one adopts the usual assumption about the noise term: (1) the 
noise is random Guassian noise which is not correlated temporally or spa- 
tially; (2) the noise is not correlated with the GBR signal, then 

Cobs = CCBR, (20) 

the three point correlation calculated from the raw observational data will 
reflect directly the three point temperature correlation of the GBR, even if 
the noise term is comparable to the signal. This is the another advantage to 
using the three point function to test Gaussianity of the initial perturbations. 
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FIGURE CAPTIONS 



Fig. 1: The dependence of the 3-pt temperature correlation on the power 
spectra index and angular separation. The solid line is for a scale invariant 
density perturbation; the dash line is for a "tilted" spectrum with spectra 
index = 0.7. 

Fig. 2: The angular dependence for the three point function in different 
scenarios: the solid line represents the 3-pt function generated by inflation; 
the dash line represents the one generated by a cosmological phase transition. 
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